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ABSTRACT 

We compute the contributions due to gluon and heavy-quark fragmentation to quarko- 
nium production at large transverse momentum in pp and pp colhsions. For inclusive 
J / if) production, there is a large contribution from the g Xc fragmentation. At large 
Pt, this is comparable to the conventional charmonium model prediction via gluon 
fusion at the ISR and dominates over the latter at the Tevatron energy. This may help 
to explain, at least partly, the large J/ ip production cross-section recently observed by 
the CDF experiment. However, the fragmentation contribution to tp' production is not 
large enough to explain the corresponding CDF data. We also present the results for 
T production at ISR and Tevatron energies. 
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The production of quarkonia in hadron-hadron and lepton-hadron collisions has 
been the subject of several theoretical and experimental investigations. Experimentally, 
the leptonic decay channels of quarkonia provide the cleanest signals of heavy quark 
production. From the theoretical point of view, quarkonium production provides a 
very important test of perturbative QCD. Both in lepton-hadron and in hadron-hadron 
collisions, the quarkonia are produced mainly through gluon-initiated subprocesses - 
i.e. photon-gluon and gluon-gluon fusion leading to the heavy-quark pair production. 
Consequently these processes are also important probes for the gluon distribution in 
the hadron. 

The colour-singlet model has been used, with reasonable success, to describe the 
leptoproduction and hadroproduction of quarkonia at large transverse momen- 
tum. In this QCD-based approach, one projects out from the full heavy-quark pair 
production amplitude, the part with the correct spin, parity and charge-conjugation 
assignments of quarkonium state. The QQ is required to be a colour singlet, 

which is achieved by the radiation of a hard gluon in the final state. The matrix- 
element so obtained is then convoluted with the wave-function at the origin, |-RoP (or 
its derivative, in the case of P-wave quarkonia). In this approach, it is possible to 
compute the production cross-sections for the different resonances in a given family 
of quarkonia: in particular, the P-wave cross-sections are computable and their con- 
tribution to the inclusive cross-section of the S'-state completely specified, once the 
branching ratios of the P-states to the S'-state are determined experimentally. The 
colour-singlet model is known to provide a reasonable description of data from the 
EMC and NMC experiments on large-p^- leptoproduction ^ of J/'?/' and on large-pr 
hadroproduction 0, ^ of J/?/) from ISR experiments. For a recent review of quarko- 
nium production, and for a discussion of the J / ip production data we refer the reader 
to Ref. §. 

The recent CDF data on large-pr J/i^ production in 1.8 TeV pp collisions at 
the Tevatron seem to indicate, however, a serious discrepancy : the theoretical 
predictions of the above charmonium model are well below the data, which suggests 
the existence of a new contribution to J/ip production at large pt- It has been pointed 
out by Braaten and Yuan ^ that a large contribution to quarkonium production at 
large Pt could come from the fragmentation of gluons and heavy quarks. Even though 
the fragmentation process is of higher order in the strong coupling constant, a^, it can 
dominate over the direct quarkonium production via fusion at large pt, where terms of 
the order of Pt/^q can easily compensate for the suppression due to the extra power 
of as- It is important to check whether the fragmentation process can explain the 
discrepancy between the theoretical predictions for J/ ip production and the data from 
CDF. 

While it is important to study the fragmentation contribution at the Tevatron 
energy in the light of the new CDF data, it is also interesting to study the magnitude 
of this contribution at lower energies. This is particularly important in view of the 
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fact that the low-energy J/ip production data has served as a benchmark for signals 
of quark-gluon plasma (QGP) formation such as J/ip suppression 0. Indeed, such 
a suppression of J/ip production has been observed [|1^ by the NA38 experiment in 
heavy-ion collisions at SPS energies (-y/i ~ 20 GeV). One important feature of the 
QGP signal is the pt dependence of the suppression; therefore, an estimate of the 
magnitude of the fragmentation contribution to J/ip production at these energies is 
important. 

In this letter, we present a complete perturbative QCD calculation of J/ip pro- 
duction at Tevatron energy, which includes the colour-singlet model prediction for 
the fusion contribution as well as the fragmentation contribution. We find that the 
gluon fragmentation contribution is large and can help at least partially to resolve the 
discrepancy between the CDF data and the earlier theoretical predictions. We also 
study these contributions for pp collisions at ^/s = 63 GeV and compare them to the 
data from the ISR experiment [O . Results for production at Tevatron and at ISR 



energies are also presented. In this case the fragmentation contribution is not large 
enough to explain the discrepancy between the fusion contribution and the CDF data 
at the Tevatron energy. This may indicate a sizable nonperturbative contribution to 
charmonium production, presumably due to the small charm quark mass. Therefore, 
we present the analogous results for T production, which may be compared with the 
future T production data from the Tevatron. 

After this work was completed, results for J/ ip production at Tevatron energy were 
presented in Refs. ||12|, |13[. Ref. [|l^ has also presented results for ip' production at 
Tevatron energy. Our results for J/ ip and ip' production at the Tevatron energy are in 
agreement with the results in Refs. [|1^, |13|. However, as emphasised above, we have 



also studied the energy dependence of the fragmentation contribution - in particular, 
we have studied J/ip and ip' production at ISR energies. We also present results for T 
production at ISR and Tevatron energies. 

In computing the inclusive J/ip px distributions, for the fusion and the fragmen- 
tation processes, the contributions of the P-wave Xc states have also to be considered. 
The x's are produced copiously and the Xi and X2 have a sizeable branching into the 
J/iIj + '~^ mode. The large-p^ production cross-section for the fusion process is given as 

E^{AB-^{J/i^,X^)X) = 
dp-^ 

[ dXiXiGa/AiXl)x2Gb/Bix2)-- ^ ^"^("^^ iJ/^,Xi)c)- (1) 

J vr 2x1 — XTcy at 

In the above expression, the sum runs over all the partons contributing to the sub- 
processes ah ( J/'?/', Xj)c; Ga/A and Gfe/^ are the distributions of the partons a and 
h in the hadrons A and B with momentum fractions Xi and X2i respectively. Energy- 
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momentum conservation determines X2 to be 



xiXtc y -2t 

X2 = — — — , (2) 

2xi — xrey 

where r = M'^/s, with M the mass of the resonance, s the centre-of-mass energy and 
y the rapidity at which the resonance is produced. 



2Mt 2pT 

xj, = Jx^T + At = xt = —^ (3) 

The expressions for the subprocess cross-sections, da /dt, are exphcitly given in Refs. 
and |]14|. We obtain da/dpx from Eq. [l| using 

^ da 1 1 da 

E = . (4) 

dp^ TT 2pT dydpT 

The fragmentation contribution is computed by factorising the cross-section for 
the process AB — ^ {J X-i)-^ i^^o ^ containing the hard-scattering cross-section 
for producing a gluon or a charm quark and a part which specifies the fragmentation 
of the gluon (or the charm quark) into the required charmonium state, i.e. 

da{AB {J/i),Xi)X) = Y, j dz da{AB cX)D^^^j/^^^^){z, fi), (5) 



where c is the fragmenting parton (either a gluon or a charm quark) and the sum in 
the above equation runs over all contributing partons. D{z,^) is the fragmentation 
function and z, as usual, is the fraction of the momentum of the parent parton carried 
by the charmonium state. The fragmentation function is computed perturbatively at 
an initial scale /io which is of the order of rric- It is then evolved to the scale typical 
of the fragmenting parton which is of the order of Pt/z, using the Altarelli-Parisi 
equation: 

dy z 

/^o-A-.(JMx.)(^) — ^*j(-'/^)^j-(JMx.)(2/)' (6) 

p j ^ y y 

where the Pij are the splitting functions of a parton j into a parton i. We consider the 
fragmentation of gluons and charm quarks alone. In principle, the contribution of the 
light quarks should also be considered but their contribution to the fragmentation at 
large pt is small and can be neglected. Further, in the evolution we consider only the 
Pgg contribution in evolving the gluon fragmentation function and the Pec contribution 
in evolving the charm quark fragmentation. The effect of the non-diagonal splitting 
function contributions can be safely neglected. The full set of initial fragmentation 
functions that we need to obtain the J/ip and the x contributions have now been 
computed. These are Dg^jj^ Dg^^ [^, Dc^^ ||l6l and -Dc^x 0' HI • 
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For the fragmentation process, the cross-section is given by a formula similar to 
Eq. I but with an extra integration over or equivalently over X2- We have 



E^iAB^{J/^,X^)X) 



J2 [ dxidx2Ga/A{xi)Gb/B{x2)Dc-,(j/^^y^^){z) — -^{ab-^cd), (7) 
with z given by 



2 X2 Xi 

For da/di(ab — > cd), we have used the lowest-order expressions. 

Using the formalism described above, we can compute the direct (fusion) and 
fragmentation contributions to J/ip production. To do so, we have to choose a set of 
parton distributions that are compatible with all the available information on structure 
functions. For J/ip production at ISR energies this is not very crucial; however, at 
Tevatron energies the values of x probed are small, and one should use a set of parton 
densities which are compatible with the low-x structure functions measured at HERA. 
In our computations, we have used |]l9l the updated MRSD-' parametrisations ||20 



for the parton densities in the nucleon, which are compatible with the most recent 
data on structure functions from HERA. The parton densities are evolved to a scale 
= /i^/4, where fi is chosen to be My for the case of direct J/ip production, and equal 
to pji'^ = pt/z for J/ip production via fragmentation. The fragmentation functions 
are evolved to the scale pr/z. We use the above choice of parameters for all the 
results presented below, except when we vary these parameters to see the sensitivity 
of our results to this choice. Another uncertainty that enters the normalisation of the 
cross-section predictions is that due to the wave-function at the origin, _Ro, and the 
derivative for the P-states, R'^. In the computation of the fusion contribution these 
appear in the subprocess cross-sections, whereas for the fragmentation contribution, 
the fragmentation functions at the initial scale are proportional to these wave-function 



factors. In Ref. |T5[ the fragmentation function for (7 — > x is written in terms of two 
parameters Hi and ifg, where Hi is related to R'l, and H'^ is a parameter that describes 
the g ^ X fragmentation via a colour-octet mechanism. The parameter H^ is rather 
poorly determined. For the parameters Rq, Hi and H'^, we have used the values quoted 
in Refs. §, |T§ (R^ = 0.8 GeV^, Hi = 15.0 MeV, H^ = 3.0 MeV). 

The results for inclusive J/ip production for both the fusion and fragmentation 
contributions at ISR {^/s = 63 GeV) and Tevatron {^/s = 1.8 TeV) have been shown 
in Fig. 1. To be able to compare directly with the data from pp collisions at the ISR, 
we have computed the invariant cross-section BEda/dp^ at ?/ = as a function of 
Pti where B is the J/ip branching ratio into leptons {B = 0.0594). For pp collisions 
at the Tevatron, we present results for Bda/dpx integrated over a pseudo-rapidity 
range |?7| < 0.5. At the Tevatron energy, the 6-quark contribution to J/ip production is 
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substantial, but this contribution has been removed from the data by the use of a micro- 
vertex detector in the CDF experiment [0. The fusion contribution (including both 
J /if) and X contributions, shown as the dashed line in Fig. 1) is significantly below the 
Tevatron data. At this energy, the gluon fragmentation contribution (again, including 
g ^ J /ip and g —>■ X shown as the dashed-dotted line in Fig. 1) is the dominant 
contribution over almost the whole range of pr values considered. The charm-quark 
fragmentation (shown as the dotted line) is insignificantly small. The sum of all three 
contributions (which is shown as the full line) is consistent with the data within a factor 
of ~ 2. The comparison with the ISR data also reveals some interesting features. The 
gluon fragmentation contribution is now smaller than half the fusion contribution at 
the lowest values of px considered, but grows steadily with increasing p^- Addition 
of the fragmentation component is seen to improve the agreement with the ISR data 
significantly. 

In Fig. 2, we show the magnitude of the variation of the total contribution, due to 
changing the scale from fi/2 to 2/i. These correspond to the upper and the lower lines 
respectively. It should be mentioned here that the former choice is favoured by several 
QCD tests. We have also studied the effect of varying the choice of parton densities 



by using GRV densities ||2T| instead of MRSD-', since both are compatible with the 
most recent structure function measurements. The effect of this change on the total 
J/ip cross-section is not significant. 

Our results for ip' are shown in Fig. 3, again for both the energies. For ip', the P- 
state contributions are, of course, absent. At the lower ISR energy, the fragmentation 
contributions are much too small so that the total contribution is equal to the direct 
contribution. For the Tevatron energy, we find that the charm quark fragmentation 
contribution dominates at large pr- The data is also contaminated at this energy 
by the 6-decay contributions. We use the 6-inclusive data from CDF and subtract the 
theoretical 6-quark contribution to get the data points shown in Fig. 3. But the sum 
of all the contributions is clearly not anywhere close to explaining the data. This is 
clearly an indication that some new mechanism is needed to explain the tp' data. It 
may be relevant to note here that the ip' has only the S'-wave contribution, whereas 
the ip is dominantly produced from the decays of the P-wave x states. A significant 
contribution from any new mechanism to S-state production will therefore lead to a 
relatively larger enhancement for ip' compared to ip. In view of the modest value of the 
charm quark mass one may speculate such a contribution to arise from nonperturbative 
effects. 

Fig. 4 shows our predictions for the case of T production. As in the case of J/ip 
production, we have included the contributions due to the decays of the P-states. 
At the ISR energy, the T comes largely from direct production via fusion while the 
fragmentation contributions are negligibly small. For the Tevatron energy, the gluon 
fragmentation contribution starts dominating only beyond a pt of 15 GeV; the rela- 
tively larger value of this cross-over point reflects the relatively larger mass of bottom 
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quark compared to charm. This also imphes that any nonperturbative contribution 
will be relatively small for the T. Thus it will be important to test this perturbative 
QCD prediction with the T production data from the Tevatron. 

To conclude, we have studied both the fusion and fragmentation contributions 
to J/^, ^' and T production at Tevatron and ISR energies. We find that at the 
Tevatron energy, the J/ip^s are produced dominantly via gluon fragmentation. The 
sum of the J/ip production cross-sections from the fusion and fragmentation processes 
is compatible with the data from CDF to within a factor of ~ 2. For i//, however, even 
the addition of the fragmentation contribution is not sufficient to restore agreement 
with the CDF data. This may indicate a significant nonperturbative contribution to 
charmonium production, presumably due to the modest value of the charm quark mass. 
Therefore it will be important to test the predicted T production cross-section with 
the future Tevatron data. 

We thank Drs. G. Altarelli and M. Mangano for several helpful discussions. One 
of us (DPR) acknowledges the hospitality of the Theory Division, CERN, during the 
course of this work. 
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Figure captions 



Fig. 1 Upper figure: The cross-section BEda/dp^ (at y = 0) for tlie process pp 

J/ipX as a function of pr at y/s = 63 GeV. The data are taken from Ref. ]TT | 



Lower figure: The cross-section Bda/dpT (integrated over the pseudorapidity 
range —0.5 < r] < 0.5) for the process pp J/ipX as a function of px at 
^/s = 1.8 TeV. The data are taken from Ref. J^. The different curves (in both 
the upper and the lower figures) correspond to the direct production via fusion 
(dashed hne), the gluon fragmentation contribution (dashed-dotted hne), the 
charm quark fragmentation term (dotted hne) and the sum of all contributions 
(solid line). 

Fig. 2 The scale dependence of the J / if) cross-section as a function of px for pp collisions 
at y/s = 63 GeV (upper figure), and for pp collisions at i/i =1.8 TeV (lower 
figure). The solid curve in both figures is for the scale /i/2 and the dashed curve 
is for the scale 2/x. 

Fig. 3 The tp' cross-section as a function of pt for pp collisions at y/s = 63 GeV (upper 
figure), and for pp collisions at y/s = 1.8 TeV (lower figure), with the different 
contributions shown as in Fig. 1. The ^/s = 1.8 TeV data are obtained as 
explained in the text. 

Fig. 4 The T cross-section as a function of pt for pp collisions at ^/s = 63 GeV (upper 
figure), and for pp collisions at ^/s = 1.8 TeV (lower figure), with the different 



contributions shown as in Fig. 1. The ^/s = 63 GeV data are taken from Ref. |jTT 
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